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Abstract 

The D0 detector is used to study pp collisions at the 1800 GeV and 630 GeV 
center-of-mass energies available at the Fermilab Tevatron. To measure jets, the 
detector uses a sampling calorimeter composed of uranium and liquid argon as 
the passive and active media respectively. Understanding the jet energy calibration 
is not only crucial for precision tests of QCD, but also for the measurement of 
particle masses and the determination of physics backgrounds associated with new 
phenomena. This paper describes the energy calibration of jets observed with the 
D0 detector at the two pp ccnter-of-mass energies in the transverse energy and 
pscudorapidity range Et > 8 GoV and |?/| < 3. 
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1 Introduction 



Jet production is the dominant process in pp collisions at -y/s=1.8 TeV. Be- 
cause almost every physics measurement at the Tevatron involves events with 

jets, an accurate energy calibration is essential. Currently, the jet energy scale 
is still the major source of systematic uncertainty in both the D0 inclusive 
jet cross section and top quark mass measurements. 

This paper describes the determination and verification of the jet energy cali- 
bration at D0. The methods developed here are based on previous work [1,2]. 
The calibration is performed for jets reconstructed with different algorithms, 
from data taken in pp collisions at = 1800 GeV and 630 GeV, but only 
representative plots are shown. For detailed information see Refs. [3,4]. 

The calorimeters arc the primary tool for jet measurements at D0. A brief 
summary of the characteristics and performance of these subdetectors is given 
in Section 2. The definition of a jet and a description of the algorithm used 
in jet reconstruction are included in Section 3. Sections 4-9 enumerate the 
different corrections involved in the jet energy calibration, explain the meth- 
ods used in their derivation, and display the numerical results. Section 10 is 
dedicated to verification tests based on a Monte Carlo derived correction. The 
results are summarized and conclusions are drawn in Section 11. 



2 D0 Calorimeters 

The D0 Uranium-Liquid Argon sampling calorimeters [5] are hermetic, uni- 
form, and provide nearly complete solid angle coverage. The Central (CC) and 
End (EC) Calorimeters contain approximately 7 and 9 interaction lengths of 
material, respectively, ensuring containment of nearly all particles except high- 
Pt muons and neutrinos. The intercryostat region (IC), between the CC and 
the EC calorimeters, is covered by an intercryostat detector (ICD) and mass- 
less gaps (MG) [5]. The ICD consists of an array of scintillator tiles located on 
the EC cryostat wall. The MG are separate single-cell structures installed in 
the CC and EC calorimeters between the module endplates and the cryostat 
wall. Given their fine segmentation and excellent energy resolution, the D0 
calorimeters are especially well suited to jet measurements and the determi- 
nation of missing energy in the plane transverse to the beams (-^t). 

A general view of the D0 calorimeters is displayed in Fig. 1, showing the 
different electromagnetic and hadronic sections. The cells are constructed of 
uranium absorber plates interleaved with Cu/GlO readout boards. The gaps 
are filled with LAr as the active medium. A cross section of a unit cell is shown 
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in Fig. 2. The coarse hadronic sections are made of copper or steel instead of 
uranium. A side view of the D0 calorimeters is shown in Fig. 3. D0 defines 
a coordinate system with the origin at the geometric center of the detector 
and the z axis along the direction of the proton beam. The transverse plane 
is defined by the cartesian axes x and y. In spherical coordinates, the three 
coordinates are described by the radius r, the azimuthal angle (/?, and the 
polar angle Q. Pseudorapidity is defined as r/ = — ln[tan(^^/2)]. The segmen- 
tation in A77 X space is 0.1x0.1, or 0.05x0.05 at shower maximum (third 
electromagnetic layer) . 

The single particle energy resolutions for electrons (e) and charged pions (tt) 
were measured from test beam data, and parameterized with a {(Je/EY = 
S"^ / E + C"^ functional form. For electrons, the sampling term S is 14.8 (15.7)% 
in the CC (EC), and the constant term C is 0.3% in both the CC and EC. For 
pions, the sampling term is 47.0 (44.6)%, and the constant term is 4.5 (3.9)% 
in the CC (EC) [6]. 

The energy of a jet measured in the calorimeters is distorted by uninstru- 
mented regions, phenomena which affect the response, energy deposits in the 
liquid argon due to uranium decays, spectator interactions, reconstruction 
and resolution effects. The calorimeter response to the different types of par- 
ticles is the most important aspect of the jet energy scale correction. Elec- 
tromagnetically interacting particles, like photons (7) and electrons, deposit 
most of their energy in the electromagnetic (EM) section of the calorimeters 
(~ 20.5 radiation lengths thick). Hadrons, by contrast, lose energy primarily 
through nuclear interactions. The showers they cause are much longer than 
those from electromagnetic particles, and extend through the three sections of 
the calorimeter: EM, fine hadronic (FH), and coarse hadronic (CH). In gen- 
eral, the calorimeter response to the EM (e) and non-EM {h) components of 
hadron showers is not the same because there are different physics processes 
involved. N on- compensating calorimeters have a response ratio e/h > 1, and 
suffer from non-Gaussian event-to-event fiuctuations in the fraction of the 
energy lost through 7r° production (/cm)- In addition, (/cm) increases logarith- 
mically with energy. As a result, such calorimeters give a non-Gaussian signal 
distribution for monoenergetic hadrons. Moreover, the signal is not linearly 
proportional to the hadron energy {non-linearity), and e/ir (ratio of electron 
to pion response) is energy dependent [7]. 

The D0 calorimeters are nearly compensating, with an c/tt ratio less than 
1.05 above 30 GeV. Figure 4 shows the e/n ratio from a Monte Carlo [8,9] 
simulation (open squares) and from test beam data (full circles) . The calorime- 
ters are also linear with energy, except at very low energy, as illustrated in 
Fig. 5 from test beam data. Figure 6 shows the Gaussian behavior of the 
variable {Eti — Et2)/{Eti + Et2) for dijet events with average jet transverse 
energy Et of 140 GeV observed in the D0 calorimeters, where Eti and Et2 
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are randomly ordered from the pair of leading transverse energy jets. This 
Gaussian signature is characteristic of the hermeticity and linearity of the D0 
calorimeters. 



3 Jet Measurement 



Jets emerging from a hard interaction deposit energy in numerous cells in 
the calorimeter. To determine the jet energy, an algorithm is needed which 
assigns calorimeter cells to jets. The final measured jet energy depends on 
this algorithm or jet definition. In this paper, only a fixed cone algorithm 
is considered [10], with jet centroid and cone size defined as {rj^et^^jet) 
n = ^J{(p- (fietY + {V- Vjetf, respectively 

The process of reconstructing jets in the calorimeter is iterative. First, towers 
[At] X Aip = 0.1 X 0.1) containing Et of 1 GeV or more are used as seeds for 
finding prcclusters, which are formed by adding neighboring towers within a 
radius of 7^ = 0.3 to seed towers. The rji of each tower is calculated with respect 
to the interaction vertex, which is determined from the tracking system [5]. 
Next, a fixed cone of radius TZ is drawn around each precluster centered at 
its centroid. A new center is then calculated using the Snowmass Accord [11] 
definitions for both the jet angles and transverse energy: 

i 

'Piet = J2^Tifi/ETiet, (2) 
i 

Viet = X EnVi/Eriet ■ (3) 



The sums are over all towers contained in the cone. This process is repeated 
iteratively, always using the Snowmass Accord definition, replacing the previ- 
ous seed or jet direction by the current one until the center becomes stable. 
In the final step, ?7jet and </7jet are recalculated using the following definition: 



(/9jet = tan 1 EyJ X E^}j 



(4) 



^jet = tan-^ L E + E EIJ E (5) 



where Ex^ — E^ sin 9i cos ipi, Ey. — E^ sin 9i sin (p^, and E'^. = Ei cos 9i. The 
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transverse energy of a jet is still the sum of the transverse energies in each 
calorimeter tower inside the cone. Those jets with Ex < 8 GeV are discarded. 
Overlapping jets are merged into a single jet if more than 50% of the of 
the lower energy jet is contained in the overlap region. Otherwise, the energy 
of each cell in the overlap region is assigned to the nearest jet. 

The particle level jet energy is defined as the energy of a jet found from final 
state particles using a similar cone algorithm to that used at the calorimeter 
level. The particle jet is never observed in data because it gets distorted when 
it interacts with the calorimeter material. But it can be clearly defined in a 
Monte Carlo sample by clustering together the final state particles, instead of 
calorimeter cells, into a fixed cone size. The particle jet includes only particles 
arising from partons participating in the hard scatter. The purpose of the jet 
energy scale calibration is to correct observed jets to the particle level. 

3.1 Jet Energy Scale 

The conversion from electronic ADC counts to energy and the layer-to-layer 
calibration are performed at the data reconstruction level and are based on test 
beam data. Detailed information on calorimeter performance and calibration 
can be found in Ref. [6]. 

This paper describes an in situ calibration based primarily on reconstructed 
collider data. This calibration relates the observed jet energy to that of the 
final state particle jet on average. The linearity and hermeticity of the D0 
calorimeters provide Gaussian response functions and, therefore, allow sepa- 
ration of the jet energy resolution from the jet energy scale correction. Thus, 
the following sections focus on the scale correction and leave aside resolu- 
tion effects. The D0 jet energy calibration is based almost entirely on data 
measurements, exploiting conservation of transverse momentum through an 
accurate determination of the event 

The particle level or true jet energy Ej^^^ is obtained from the measured jet 
energy E"^^^ using the following relation: 

pptci ^ E^--Eo{n,r^,C) 

it:jet(7^,77,£;)5(7^,77,£;)' 

where: 

• Eo{TZ,r], C) is an offset, which includes the physics underlying event, noise 
from the radioactive decay of the uranium absorber, the effect of previous 
crossings (pile-up), and the contribution of additional pp interactions. The 
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physics underlying event is defined as the energy contributed by spectators 
to the hard parton interaction which resuhed in the high-p^- event. The 
offset increases as a function of the algorithm cone size TZ. The dependence 
on luminosity, £, arises from the contribution of additional interactions. Eq 
is also parameterized as a function of rj, measured with respect to the vertex 
of the interaction (physics 77). 

• i?jct('7^, i]j is the energy response of the calorimeter to jets. It is parame- 
terized as a function of jet energy after the offset subtraction. Rjet is inde- 
pendent of the algorithm cone size; however, its parameterization versus jet 
energy is cone size dependent because larger cones encompass a larger frac- 
tion of the energy of the cluster. Because the various detector components 
are not identical, i?jct is also dependent on detector pseudorapidity, that is rj 
calculated from the geometric center of the detector and not the interaction 
vertex. Rj^i is typically less than one, due to energy lost in uninstrumented 
regions between modules, differences between electromagnetic and nuclear 
interacting particles {e/h > 1), and module-to-module inhomogeneities. 

• S{TZ, 7], E) is the fraction of the particle jet energy that is deposited inside 
the algorithm cone. The jet energy is corrected back to the particle level, 
and therefore calorimeter showering effects must be removed. S is less than 
one, meaning that the effect of showering is a net flux of energy from the 
inside to outside the cone. S depends strongly on the cone size 7^, energy, 
and 7]. It is parameterized as a function of jet energy after offset subtraction 
and response correction, and binned in terms of physics rj. 

Most of the examples presented in this paper correspond to jets recon- 
structed with SkTZ— 0.7 fixed cone algorithm. Jets reconstructed with other 
cone sizes or at ^/s = 630 GeV were also calibrated. This information is 
included in Refs. [3,4]. 



4 Offset Correction 

The pp inelastic collisions are classified as non-diffractive, single diffractive 
and double diffractive. An inelastic non-diffractive collision, also sometimes 
called hard core interaction, is an event in which both the proton and the 
anti-proton break up. A particular case of a hard core interaction is high-p^ 
or hard parton scattering. 

The offset correction is determined for non-diffractive events. It subtracts en- 
ergy which is not associated with the high-p^- interaction. This excess energy 
includes the contributions of uranium noise, pile-up, and additional pp interac- 
tions. Pile-up is the residual energy from previous pp crossings, which results 
from the long shaping times of electronic signals in calorimeter readout cells. 
Pile-up depends on the number of interactions in the previous beam cross- 
ings, and is therefore luminosity dependent. The underlying event includes 
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soft interactions between spectator partons that constituted the colUding pro- 
ton and anti-proton. This energy can be subtracted or not, depending on the 
characteristics of each particular physics measurement. 

A crossing selected by a high-pr trigger can be modeled as the sum of a hard 
parton scattering overlaid with a zero bias (ZB) event at the same luminosity. 
The ZB trigger accepts every pp crossing, regardless of whether or not it 
contains an actual collision. The average number of hard interactions in a 
high-pr event, (A^), may be written as: 

(AT) = 1 + (ATzb) . (7) 

Here, (A^zb) is the average number of hard core interactions in a ZB event and 
can be expressed as: 

oo 

(7V,3) = ^nP(n), (8) 

n=0 

where the probability of having n hard core collisions at a given luminosity, 
P(n), follows a Poisson distribution. 

The total offset correction, Eq, can be presented as the sum: 

noise + -Spile 

= Eue + (A^ZB)-E'ue + -E'noise + -E^pile , (9) 

where E^^ is the energy associated with the physics underlying event, {NzB)Eue 
is the contribution to the underlying event from additional pp interactions, 
E'noise is the energy from uranium decay, and EpHe is the pile-up energy. For 
the actual correction, the last three terms may be combined, yielding: 

^O = ^ue + ^0- (10) 

4-1 Physics Underlying Event 

The energy from spectator partons is measured using low luminosity minimum 
bias (MB) data {C = 0.1 x 10^° cm~^sec~^). The MB trigger requires an 
inelastic collision in which both the proton and the anti-proton break up. MB 
events are, therefore, dominated by hard core interactions. The E^ density 
per increment At] x A(p in MB events, i^MB, is a measure of the transverse 
E^e + -Snoise + ^'piie density. The contribution of uranium noise and pile-up. 
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-^ZB ) is measured in a low luminosity ZB sample with the requirement that 
there is no hard interaction. D,,e is then: 



It depends both on the center-of-mass energy and pseudorapidity. The energy 
to be subtracted from a jet, E-^^, is D^^ x A^^^, where Arj_^ is the area of the 
jet in Tj-ip space. Figure 7 shows the dependence of D^^ on physics r] for both 
^/s — 1800 GeV and ^/s — 630 GeV coUisions. The uneven shape results from 
imperfect calibration of the intercryostat detector and loss of EM coverage 
in the intercryostat region. Because the physics underlying event is associated 
with the soft interactions in a single pp collision, it is independent of luminosity 
and of the number of pp interactions in the event. 



4-2 Uranium Noise, Pile-up, and Extra Interactions 

The contribution to the offset from uranium noise, pile-up, and additional pp 
interactions, Dq, is the measured Et density in a ZB sample. Figure 8 displays 
parameterizations of Dq as a function of physics r) for different luminosities. 

Dq also depends on the occupancy, defined as the number of calorimeter cells 
which were read out after zero-suppression divided by the total number of 
cells in the jet. Zero-suppression is the mechanism that suppresses cells with 
energy depositions inside a 2(T window around the average pedestal value. 
Therefore, Dq from the ZB data must be extrapolated to a value consistent 
with average occupancies for jet data. This difference in occupancy dominates 
the uncertainty in the offset correction. It contributes an error of ± 0.25 GeV 
to Dq. In addition, statistical and systematic errors in fitting the data are 
approximately 8%. 



5 Response: The Missing Et Projection Fraction Method 

The jet reconstruction algorithm maps charge collected in the liquid argon to 
energy. It is based on single particle test beam data, assuming ideal instru- 
mentation and a linear response. The overall response after reconstruction is, 
however, less than unity, due to a non-linear response to low energy particles, 
dead material, and module-to-module response fluctuations. This section de- 
scribes the missing Et projection fraction method [13] developed to measure 
the calorimeter response to jets in reconstructed data. 
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5. 1 Definitions 



D0 makes a direct measurement of the jet energy response, using conservation 
of pt in photon-jet (7-jet) events. These events are composed of one photon 
balanced in pt by one or more jets. The missing transverse energy, is 
defined as: 

$T = -{^E,.,Y.Ey)l , (12) 



where Ex^^ Ey^ are defined as in Section 3 and i runs over all calorimeter 
cells with |?7| < 4.5. This definition implies that the energy deposit in each 
cell may be treated as though it were deposited by a zero mass particle, so 
that $j,='^j,. In an ideal calorimeter, any non-zero indicates the presence of 
particles which failed to deposit their energy in the detector, such as neutrinos 
or high-pT muons. For a photon-jet event in a real detector, a non-zero 
measures the overall imbalance of transverse energy in the calorimeter due 
to differences in response to photons and jets. This can be used to measure 
the calorimeter response to jets, i?jet, relative to the precisely known photon 
response. 

In 7-jet events, the particle level or true photon and recoil transverse energies 
Et^ and i^rrecoii, satisfy: 

Er-y + -ETrecoil = . (13) 



In a real calorimeter, however, the photon and jet responses (i?em and i?recoii) 
are both less than unity, and the equation is modified to: 

-, -t -> meas 

iTimeas i zrimeas jp ( -\ A\ 

Ej, + -C/Trecoil = "-^T > l^^j 



where E^^^ = RcmET^ and i?™*ccoii = -Rrccoii-E'Trccoii- The energy scale for elec- 
tromagnetically interacting particles is determined [14] from the Z e~^e~, 
J/ip, and 7r° data samples, using the known masses of these resonances. If 

— * 

^meas jg corrected for energy scale in the 7-jet data sample, Eq. 14 transforms 
into: 



Et'y + -Rrecoil-E'Trecoil = ~Pt 
Et^ -\- Rrecoil'f^T-y ' -^rrccoil = ^^T^ ■ 

1 , D ■ -^Trecoil flxj ■ $t n k\ 

i + -Krecoil — ^ , [i-O) 
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where 77,^7 = ET'y/\ET'y\, and $j, is the missing transverse energy recalculated 

— * — * 

after the photon correction. Eq. 13 can be written as Ex-y — —nr-y ■ -E'Trecoib 
then -Rrecoii is: 

i?.ecoiI = 1 + %^ . (16) 

-C/T7 



In the special case of a 7-jet two body process, and in the absence of offset 
and showering effects, i?rccoii would be the ratio E^^-'^f / E^^^^ of measured to 
particle jet transverse energies. In the presence of offset and showering losses, 
-Rrecoii IS the energy response of the calorimeter to jets, Rjet, where "jet" is 
the leading jet of the event. This is a good approximation if the difference in 
azimuth between the 7 and the leading jet is close to tt. 



5.2 The Energy Estimator, E' 



Rjet is measured at D0 as -Rrecoii from a 7-jet sample, using conservation of 
transverse momentum (the terms transverse momentum and transverse energy 
are used interchangeably). The response is, however, dependent on jet energy 
rather than its transverse component. This is because c/tt and the particle 
composition of jets are energy dependent. 

Measuring Rjet directly as a function of E^^^ is problematic. The finite jet 
and photon energy resolutions, steeply falling photon cross section, trigger 
and reconstruction thresholds, and event topology contribute biases which 
must be removed. Most of these biases and smearing effects are reduced to 
negligible levels by binning the response not in EJ^^^ but instead in some 
better measured quantity which is strongly correlated with E^^^^^. For this 
quantity the jet energy estimator E' is chosen and defined as: 

E' = Et^ ■ cosh(7^) , (17) 



where Etj includes the electromagnetic scale correction. Both Etj and the 
jet pseudorapidity are measured with high resolution compared to EJ^^^ . 
After binning the response in terms of E', the dependence of Rjet on EJ^^'^^ is 
obtained by measuring the average E^^'^^ in each E' bin. Figure 9 illustrates 
the E' —>■ E^^^ mapping procedure. 
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5.3 Tests of the Method 



The method is tested using a parametric simulation. The simulation gener- 
ates 7-jet events according to a cross section with a given Et dependence. 
The energies of the photon and the jet are smeared and scaled with energy 
resolutions and responses measured from data. 

Figure 10 (top) shows a comparison between Rj^t determined from the sim- 
ulation directly as a function of E^^^ (open circles), and the same quantity 
binned in terms of E' and then mapped onto E"^^'^^ (closed circles). While the 
closed circles are in very good agreement with the input response (solid line), 
the open circles exhibit the smearing effect arising from a finite jet energy 
resolution. Figure 10 (bottom) shows excellent agreement between the fits to 
Rjet versus E' and the input function given different reasonable assumptions 
for the Et dependence of the photon cross section, such as E^^,E:^'^, and E^'^. 



6 Response: Photon Sample Selection 

A large event sample is necessary to measure the jet response of each calorime- 
ter with a reasonably small uncertainty over a large range of energies. Thus, 
Rjet is measured from a collider data sample obtained by applying a set of cuts 
which retain a large sample while limiting the systematic biases. For the pur- 
pose of cahbration, the "7-jet" sample is not restricted to direct photon events. 
It also includes events with one jet well contained in the EM calorimeter and 
little nearby hadronic energy. 

6.1 General Cuts 

From the events with photon candidates accepted during data collection, only 
those with one or more reconstructed EM clusters are used. Events with noisy 
cells which fake or distort jets are rejected. The Main Ring, the injector for the 
Tevatron, passes through the hadronic calorimeter and is active during data 
taking for p production. Events containing particles related to Main Ring 
activity are, therefore, removed from the sample. The muon system is used 
as a loose veto to avoid events with bremsstrahlung photons from cosmic 
ray muons. If the leading photon Et is less than 30 GeV, there must be no 
reconstructed muons in the event; otherwise the event is accepted if muons 
have pTfj, < 100 GeV/c. Because events with primary vertices far from the 
center of the detector can distort Et measurements, the 2;-coordinate of the 
vertex must satisfy \z\ < 70 cm. 
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Fiducial cuts on detector rj are applied to ensure good containment of the EM 
cluster inside the calorimeters {\ri\ < 1 or 1.6 < \r]\ < 2.5). Candidates near 
detector boundaries are also rejected. There must be at least one reconstructed 
jet in the event, and the leading jet must be contained in the CC {\r]\ < 0.7) 
or the EC (1.8 < \r]\ < 2.5) for the CC and EC calibrations, respectively. 

There are two types of background to the calibration sample, removed as 
described below. 

6.1.1 Instrumental Background 

The instrumental background is contributed by events in which the photon, 
or pair of photons from a highly EM jet (mostly 7r° jets), is not well isolated 
from a significant amount of hadronic energy. The following variables are used 
in the selection procedure to remove background: 

• The fraction of the cluster Ex deposited in the EM layers of the calorimeters. 



• The cluster isolation in the plane transverse to the beam, /iso, defined as: 



where E'tot and Eem are the total energy and the EM energy within a cone 
radius 71. 

• The total charge in the transition radiation detector [5], St, measured in 
units of the charge of the electron. 

• The presence of a track in a road projected from the EM calorimeter cluster 
through the tracking chambers to the event vertex. Track match significance, 
(Ttrk, is defined as: 



where Aip and Az are the distances between the track and the cluster along 
the azimuthal and beam directions. 
• The ionization in the central tracking detectors, dE/ dx, in units of minimum 
ionizing particles. 

Table 1 shows the cuts on photon candidates to reject background. The three 
cuts based on tracking information, e^, crtrk, and dE/dx must be satisfied for 
events with Et^j > 25 GeV only if ^j./Et^ > 0.65. Instead, only one of the 
three is required if Er-y < 25 GeV or ^j^/Et-^ < 0.65. 

There is a bias in i?jet due to the presence of highly electromagnetic jets in 
the sample, as compared to photons. This bias is studied using a test [12], 



EMF. 



iso — 



£;tot(7^ ^ 0.4) - ^em(7^ = 0.2) 

E,^{n = 0.2) 
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which measures the consistency of the shower shape of the cahbration sample 
with that expected from an electromagnetic shower. A matrix is constructed 
with transverse and longitudinal shower shape information from the photon 
candidate, and then compared with the same type of matrix built from electro- 
magnetic clusters in a Monte Carlo sample. A cut on ^ removes a significant 
number of events but is a good discriminator between real photons and highly 
electromagnetic jets. Figure 11 shows the magnitude of the instrumental back- 
ground bias on the jet response. As the cut on tightens, the fraction of real 
photons in the sample increases, but the response changes only slightly. So 
that a large sample can be retained, a cut is not applied in the event selec- 
tion procedure, and a bias on the order of 0.7% towards a larger response is 
present. A large fraction of this bias cancels against the topology bias, as will 
be discussed in Section 8. 

6.1.2 Physics Backgrounds 

The physics backgrounds include Drell-Yan, Z — > e+e~, and W ^ eu events. 
In events with two electromagnetic clusters, such as Z ^ e'^e~ production, 
i?jet is approximately unity. The effect of the Z boson background can be de- 
termined by applying a cut on the dielectron invariant mass, rUee- This back- 
ground is very small and, therefore, the cut 75 GeV/c^ < mge < 105 GeV/c^ 
produces a negligible change in the measured response. 

W ^ eu events have a large in the direction of the neutrino, so they are 
not useful for calibration purposes. They are removed using the tracking cuts 
listed in Table 1, combined with a cut, as described before. After the cuts, 
the remaining bias on i?jet due to physics background, determined from Monte 
Carlo and data samples, is <j 0.5%. 

6.2 Topology Cuts 

The relation Rjet = -Rrecoib where the subscript jet refers to the leading jet in 
the 7-jet event, is exact only in the case of a two body process. Photon events 
may contain more than one reconstructed jet along with a number of energy 
clusters which are not reconstructed as jets. Variations in topology, therefore, 
contribute a systematic error to the response measurement. 

To remove this topology bias, the photon and the leading jet must have a 
difference greater than 2.8 radians in azimuth. The residual error is investi- 
gated by studying the change in response for different energy bins as the Acp 
cut is varied from 2.1 to 3.1 radians. Figure 12 shows the change in R^et as 
the Aip cut is tightened for the 50 < E' < 60 GeV bin. The circles corre- 
spond to an integrated distribution (lower A(p bins include the events of the 
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higher A(/? bins), while the open squares represent the differential distribu- 
tion. Requiring Atp > 2.8 radians is a good compromise l:)etween a large event 
sample and bias minimization. The residual error after the cut is estimated 
with the parameterizations to the integrated (linear fit) and differential dis- 
tributions (exponential fit). The fractional difference between i?jet at A(p — tt 
and A(p > 2.8 radians is an estimate of the remaining topology bias. This 
difference as a function of E' is plotted in Fig. 13 for both the CC and the 
EC. The plots show that, after requiring Aip > 2.8 radians, the remaining 
bias is approximately 0.5-1% towards a lower response and consistent with a 
constant dependence on E'. A large fraction of this bias cancels against the 
effect of instrumental background, as will be discussed in Section 8. 

6.3 Multiple Interaction Cuts 

Events which consist of several overlaid interactions are called multiple inter- 
action events. Usually, only the "primary" interaction will produce high-E'^' 
objects. The additional interactions can, however, degrade the accuracy of 
the vertex determination for the primary interaction. On average, jets recon- 
structed with an incorrect vertex are assigned a higher pseudorapidity, yielding 
a larger event E' and a lower jet Et- The overall effect is an increase of the 
in the direction of the jet, thus lowering the measured response. This ef- 
fect grows as luminosity increases, as a result of the increasing proportion of 
multiple interaction crossings. 

To reduce this bias, a low luminosity sample is used. This sample, however, 
still contains a number of multiple interaction events. These are removed by 
requiring the trigger, tracking detectors, and calorimeter information to be 

consistent with a single interaction. Figure 14 shows the change in the mea- 
sured response from a low luminosity sample when the single interaction re- 
quirement is applied. The low tails in the i?jet distributions are removed by 
this cut. 

To measure any residual luminosity dependence of the response after the single 
interaction cut, the response is measured as a function of luminosity for various 
E' bins for CC and EC jets, as shown in Figs. 15 and 16. A linear fit of R^et 
versus JC allows an extrapolation to zero luminosity. In Fig. 15, the CC data 
show very little sensitivity to luminosity. The EC data, however, show a slight 
trend to lower response with increasing luminosity. The estimated difference in 
response between a luminosity of 5 x 10^° cm~^sec~^ and zero is approximately 
0.5% for the EC jets. Because the sample used in this analysis is integrated 
over this range of luminosities, any residual luminosity effects are less than 
0.25%. 
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7 Response: Rapidity Dependence 

After test beam calibration, the D0 calorimeter system still shows non-uniform 
regions in pseudorapidity Given that most physics measurements need a high 
level of accuracy at all rapidities, ?7-dependent corrections become necessary. 
These arc the cryostat factor and the IC corrections applied in consecutive 
order to both the jet energy and the event 

The 77-dependent corrections are applied as multiplicative factors to the jet 
energy. Given that the energy corrections are not derived for calorimeter cells 
but for physics objects, it is not possible to recalculate jj^j. using the definition 

-± corr 

in Section 5. Instead, the corrected missing transverse energy, , is deter- 

mined by adding to the measured or reconstructed missing transverse energy, 
, the difference between the uncorrected and the corrected vector sums 
of the Et of the individual physics objects (jets and photons in the calibration 
sample). In other words: 

= + E - E--) + E (%r - • (18) 

7 jet 

This formula depends on the jet algorithm because it is based on reconstructed 
physics objects instead of calorimeter cells. It is, however, a good approxima- 
tion to the corrected for the photon scale and the rapidity dependence of 
jet response, if the jet cone size is large (7^=0.7). 

7.1 Cryostat Factor Correction 

The cryostat factor, Fcry, is defined as the ratio Rf^^/Rf^^. Fcry is expected to 
be independent of E' because the construction of the GG and the EG is similar. 

The measured response versus E' before corrections is shown in Fig. 17. -Rjet 
is plotted in three detector rj regions: GG (|?7| < 0.7), IG (0.7 < |?7| < 1.8), 
and EG (1.8 < \ri\ < 2.5). 

Fcry is measured where the GG and the EG data overlap, in the range 60 GeV < 
E' < 180 GcV. In the EG data, an additional cut Et-^/ > 25 GcV is necessary 
to remove the reconstruction threshold bias due to jet finite energy resolution. 
Fcry is measured separately in the EG north (—2.5 < 77 < —1.8) and the EG 
south (1.8 <ri < 2.5) calorimeters. The ratio F^^yZ-^cry = 0.997 ± 0.003 shows 
that Fcry is the same in both EG's within errors. Figure 18 shows the measured 
value for the cryostat factor, Fcry = 0.980 ± 0.007 (stat), and illustrates its 
derivation from a fit of a constant function to the measured RfJ^/R'^^f versus 
E'. 



20 



The most important conscqTicncc of the independence of Fcry on E' is the 
possibihty of using EC data to extend the range of the CC response measure- 
ment. This is possible because forward jets have higher energies than central 
jets with the same Et- Thus, the energy range of the response measurement 
can be extended from ~ 150 GeV to nearly 300 GeV. 

7.2 IC Correction 

The intercryostat region, covering the pseudorapidity range 0.8 < \rj\ < 1.6, is 
the least well instrumented region of the calorimeter system. The IC is a non- 
uniform region covered by different types of detectors. A substantial amount of 
energy is lost in the cryostat walls, module endplates, and support structures. 
Additionally, in the range 1.2 < |?7| < 1.4, the system lacks electromagnetic 
calorimetry and the total depth falls slightly below 6 interaction lengths. As 
a consequence, the response has a residual r] dependence in this region, even 
after the cryostat factor has been applied to EC jets. 

7.2.1 Event Selection and Method 

The IC correction to the response is measured using both 7-jet and jet-jet 
events. Because the dependence of the response on pseudorapidity is related 
to detector inhomogeneities, -Rjct is measured as a function of detector pseu- 
dorapidity. For the 7-jet sample, the photon is kept central {\ri\ < 0.5), while 
the jet Tj is unconstrained. For the jet-jet sample, one of the two leading Et 
jets is required to be central (|?7| < 0.5), and there is no restriction on the rj 
of the other jet. In this case, Rjet is defined as: 

^central j et i n N 

-rtjet = i + , (lyj 

-C/T central jet 

and is measured as a function of the jet detector rj. The 7-jet events are 
used only at low jet Et, as the sample is statistically limited at high Et- The 
correction for high Et jets is determined from the jet-jet data. For two identical 
Et jets, one central and one far forward, it is more likely for the central r) jet to 
fluctuate to a higher measured Et than the forward jet. This is because, for the 
same Et, forward jets have higher energies and, therefore better resolutions 
than central jets. In addition, forward jets are kinematically suppressed. This 
causes a bias in the jet-jet sample above I77I = 2 which is removed with a cut 
on the product of the transverse energies of the two leading jets. 

The IC correction is performed before the energy dependent response correc- 
tion. Because the energy dependence of Rjet is folded into it!jet versus 77, this 
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function is not a constant even for a uniform detector (ideal rj dependence); 
it must be, however, smooth. After the cryostat factor correction is apphed, 
the IC correction is determined by means of a smooth interpolation through 
the IC of a fit to the measured R^et versus 77 in the CC and the EC. 

7.2.2 Results 

To determine the ideal rj dependence of the response in the IC, the response 
is fit to the function Rj^t = a + b ■ ln[cosh(r7)] in the CC and EC regions. This 
functional form is derived from the energy dependence of the response, which 
is well described by: 

Rj^i^ a + b-\nE . (20) 

Thus: 

Rjet = a + b ■ \n{ET ■ cosh rj) 

^a + b-\nET + b-\n{coshr]) . (21) 

Reordering terms, and for a fixed Et bin: 

Rjet^a + b- ln(cosh 77) . (22) 

The data are fit in the range \ri\ < 0.5 and 2 < \rj\ < 2.5. A cut on Eti ■ -E't2, 
which requires this product to be greater than the Et threshold for which a 
particular jet trigger is fully efficient, is applied to avoid the resolution bias 
in jet-jet data above rj = 2. The two halves of the detector, r] > and f] < 0, 
are treated separately. Figures 19 and 20 show the measured r] dependence 
for jet-jet and 7- jet data respectively, along with the fit to an ideally uniform 
detector. 

The correction is determined in bins of 0.1 in detector rj. The difference be- 
tween the fit to the ideal response and the measured value is used to obtain a 
correction factor F^. Figure 21 shows as a function of the average central 
jet Et {TZ = 0.7) for several rj bins. Also shown is a linear fit to the data, 
which is used to determine the correction factor as a function of jet Et- 

7.2.3 Error Analysis and Verification Studies 

The accuracy of the IC correction is measured using the fact that the response 
versus 77 should be constant after all corrections are applied. This is shown in 
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Fig. 22, where the corrections apphed to i?jet in 7-jct data arc the offset, the 
cryostat factor, the IC, and the energy dependence (sec Section 8). The hor- 
izontal line indicates the ideal rj dependence of the corrected i^jet- Figure 22 
also shows that the measured i?jet is shghtly less than unity. The reason is 
that the correction to the is only apphed on found objects. As the unclus- 
tered energy is not corrected for response effects, the average total corrected 
transverse energy of the event is not expected to be exactly zero. 

Figure 23 shows the response versus 77 measured from jet-jet data after all 
corrections have been applied. The error on is determined from the residuals 
of the corrected 7-jet and jet-jet data with respect to the ideal rj dependence. 
Table 2 shows the means and RMS values for the residual distributions in 
7-jet data, and Table 3 displays the means and RMS values for the residual 
distributions derived from jet-jet data. The values measured from each sample 
are in good agreement. 

There is some evidence of a residual 1] dependence beyond l^^l = 2.5 in the 
7-jet measurements. This bias towards a lower response at high \ri\ may be 
due to the breakdown of the constant cryostat factor approximation at very 
high pseudorapidity, or to energy lost down the beam pipe. To account for 
this effect, an additional r^-dependent error is assigned. It starts from zero at 
\r]\ = 2.5 and increases linearly up to 3% at \r]\ = 3. 

The IC correction increases for decreasing cone sizes. It is determined sepa- 
rately for each cone size. Once the 77-dependent corrections arc applied to both 
the jet energy and the event ^j,^ the response becomes constant as a function 
of pseudorapidity to within ~ 1% (|7;| < 2.5). 



8 Response: Energy Dependence 

In Section 2, the ratio e/vr was shown to be greater than unity and energy 
dependent. R^ct must, therefore, also be energy dependent. The energy de- 
pendence of i?jet at low to moderate E' {E' < 100 GeV) is determined from 
low-£'r photons and CC jets (I77I < 0.7). For high E' (> 100 GeV), i?jet is mea- 
sured from EC (1.8 < \ri\ <2.5)jets, after F^y and corrections are applied, 
taking advantage of the fact that the CC and EC are similarly constructed. 
In other words, the EC measurement is normalized to the CC measurement 
using Fcry, in order to extend the energy reach of the global -Rjct measurement. 
This Section starts with a discussion of the low-F^ bias arising from recon- 
struction and resolution effects, and follows with a detailed description of the 
measurement of it!jet versus energy. 
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8.1 The Low-Et Bias 



The combination of reconstruction inefficiencies, a minimum jet Et cut, and 
tfie finite jet energy resolution produces a bias in tfie measured response. 
At the Et threshold for jet reconstruction E™™ (set at 8 GcV for all jet 
cone algorithms at D0), the jet fractional Et resolution is approximately 
50%. Thus, the net migration of \ow-Et jets to higher values is very large. 
At the same time, jets which fluctuate below are not reconstructed. 

This low-end truncation biases the average jet Et to higher values, biases 
the ]pT lower values, and therefore biases the response to higher values. 
The jet reconstruction efficiency is less than unity due to inefficiencies in the 
algorithm for turning seed towers into jets. Consequently, the \ow-Et bias is 
present above the minimum jet Et reconstruction threshold (up to 20 GeV), 
and is both cone size and r^-dependent. 

Because the jet energy response is determined from the event $ti ^T-y, and 
n-y, it may be measured without requiring any reconstructed jet in the event. 
To determine the effect of the bias, the response is measured as a function 
of i?T7 for those events both with and without at least one reconstructed jet. 
The low-£^T bias is then the ratio: 

p _ -^jel(>ijel.) 

-Kbias — • y^o) 

-'ijet(no jet required) 

An independent estimate of the 1ow-£'t bias is obtained from simulated 7-jet 
events, constructed as in Section 5. The jet response is then "measured" using 
the missing Et projection fraction method. The low-ii^T bias is determined 
from the ratio of the response from a sample where the E"™™ cut is modeled 
to the response from a sample with no such restriction. This bias depends 
on the jet energy resolution, Et™, and the reconstruction efficiency, while the 
offset and the Et dependence of the photon cross section contribute very little. 
Figure 24 compares the simulated bias and the data measurement for 0.7 cone 
jets. The results agree within the errors of the simulation (dotted band). The 
Ruas correction is obtained from parameterizations of the data measurements 
for different cone sizes, and the uncertainty is determined from the simulation. 

8.2 Response versus E' 

After application of the Iow-Et bias, offset, cryostat factor, and IC corrections, 
-Rjct is recalculated as a function of E' (Fig. 25). The average response in each 
E' bin is determined by fitting the response distribution with either a single 
or double Gaussian parameterization. A single Gaussian is sufficient for E' 
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bins with less than 100 entries. These fits provide a response measurement 
less affected hj misvertexing than the arithmetic means of the distributions 
because they are less sensitive to an occasional outlying point far from the 
peak. In Fig. 25, the points are plotted at the average value of E' within 
each bin. While the 77 dependence of Rjet is dramatically reduced, as shown in 
Fig. 25, the IC is slightly lower compared to the CC and EC regions. Although 
the CC, IC, and EC responses are the same within the error of the ry-dependent 
correction, only the most accurate CC and EC data are used to derive the jet 
response. 

After determining the response versus E', a mapping between E' and the 
average jet energy is obtained separately for each jet algorithm. The mapping 
for 0.7 cone jets is shown in Fig. 26. The observed difference in this mapping 
between EC and CC jets is consistent with showering effects, as will be shown 
in Section 9. 



8.3 Constraining High Energy Response with Monte Carlo 

The jet energy response can be determined from the data for jets with energy 
up to ~ 300 GeV, as shown in Fig. 26. However, energy calibration is required 
for the highest energy jets in the data, nearly 600 GeV. The inclusion of 
Monte Carlo information at high energy is necessary in order to constrain the 
jet response extrapolated from the available data. 

A set of 7-jet events was generated using HERWIG [15], processed through 
the GEANT-SHOWERLIB [16] fast detector simulation, and reconstructed with 
the standard photon and jet algorithms. SHOWERLIB is a library that con- 
tains single particle calorimeter showers obtained using the GEANT full detec- 
tor simulation. Because GEANT is a time consuming computer package, the 
SHOWERLIB approach is convenient. The jet response as measured from this 
sample is inconsistent with the response as measured from data, as shown in 
Fig. 27. This is because the Monte Carlo does not correctly model the e/vr 
ratio (see Fig. 4 in Section 2). Therefore, a different approach is followed to 
generate the Monte Carlo points at high Ex- 

The method consists of returning to the particle level jets in the Monte Carlo. 
For a given particle jet, the sum of the energies of all final-state particles 
in the jet cone gives the total energy of the jet. The detector may then be 
simulated by convoluting the particle jet with single particle test beam data. 
Each individual particle is scaled by the appropriate response measured from 
the test beam, parameterized in terms of energy. The response functions used 
are shown in Figs. 28a-c. Electrons and photons are scaled by the electron 
response, and 7r° mesons are scaled as two photons, each with half of the 
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energy of the pion. All other particles are scaled by the response of charged 
pions. Note that there are no test beam data below 2 GeV. Figures 28a-c show 
three different assumptions for the behavior of the response in that low energy 
range. For comparison, Fig. 28d shows the single particle responses given by 
the SHOWERLIB Monte Carlo. 

As a test, illustrated in Fig. 29, the SHOWERLIB Monte Carlo response as 
determined using the $j. projection fraction method is compared with the 
response derived from the single particle convolution method, using the Monte 
Carlo response functions (Fig. 28d). The two responses agree to within 1%. 
The shapes of the response curves derived from the three different assumptions 
shown in Figs. 28a-c are then compared with the shape of the response derived 
from data. As shown in Fig. 30, the best agreement is found from the model 
in Fig. 28b (extrapolated electron response, fixed pion response). 

The model-dependent error in the extrapolation depends upon the lower cutoff 
energy for normalizing to the data. The two extreme single particle response 
models are used: the extrapolated and the fixed responses (Figs 28a and 28c 
respectively). The Monte Carlo curves are normalized to the entire energy 
range where data exist, as shown in Fig. 31. The assumption is that these 
two extreme models bound the true response in the region where no data are 
available. The response is estimated as the mean value of the two extreme 
models at 500 GeV, and the error due to choice of model (inner error bar) as 
the difference divided by \/l2 (~ 0.7%). To this, an uncertainty of ~ 0.6% 
is added hnearly, based on the closure test illustrated in Fig. 29. The small 
normalization and fit errors are also included in the full error bar shown. 



8.4 Response Fits 

The response versus energy for 0.7 cone jets is shown in Fig. 32. The data are 
fit with the functional form: 

R^^t{E) = a + b-hiE + c-{hiEf . (24) 



A logarithmic functional form is motivated by the fact that the EM content 
of a hadronic shower slowly increases with increasing shower energy [17]. In 

the development of such a shower some fraction (/cm) of the energy is spent 
in the production of 7r° and t] mesons. A hadronic shower, therefore, has both 
an electromagnetic and a non-electromagnetic component. As defined in Sec- 
tion 2, e and h are the responses of the calorimeter to the EM and non-EM 



26 



components of a hadronic shower. If tt is the response to charged pions, then: 



where, in general, (/em) is reasonably well described by ~ a • InE' (GeV). In 
a completely compensating calorimeter e/h—1, and the ratio e/ir is unity for 
all energies, li e/h ^ 1, the ratio e/vr tends to one at high energies but grows 
exponentially as energy decreases. If e = 1, then tt would be a + b - InE. 

The logarithmic behavior of Rjet versus energy is suggested by the test beam 
measured e/ir ratio as a function of beam energy [6], and verified directly from 
the measurement of i?jet using collider data. The inclusion of the (In term 
in Eq. 24 improves the agreement of the fit with the data in the quickly varying 
low energy region, while preventing the slope of the function from rising faster 
than the data at high energies. 

The three sets of points shown in Fig. 32 correspond to CC jet data (open 
circles), EC jet data (filled circles), and Monte Carlo data (star). The three 
lowest energy points have large errors (nearly fully correlated point-to-point 
in energy), as a result of the \ow-Et bias correction. They are not used in the 
fit which includes points above 30 GeV. The of the fit is 10 for 15 degrees 
of freedom. Table 4 lists the fit parameters for the various cone algorithms. 
The calorimeter response is the same for all cone sizes as expected, although 
the parametrizations are slightly different because each algorithm associates 
a different energy to the same cluster. 



8.4- i Errors from Response Fits 

The errors in the fit parameters, hsted in Table 4, represent one standard devi- 
ation uncertainties, calculated from the — Xmin + 1 surface in the parameter 
space. However, the probability for all fit parameters to simultaneously take 
on values within the one standard deviation region (x^ = xLia + 1) is signif- 
icantly less than 68% [18]. To account for this reduced probabihty, a surface 
is mapped out in parameter space containing a region with a 68% probability 
for parameter fiuctuations in all the response fits. For three parameters this 
corresponds to the surface defined by = Xmin + 3-5 [18]. The points lying on 
this surface are then mapped back onto the response versus energy plane. At 
each energy, the error is determined by the parameter set giving the greatest 
deviation from the nominal response. The result for the 0.7 cone algorithm is 
shown in Fig. 32. 

The band represents the error on the fit to the response, and is used in place 
of one derived from the errors in Table 4. At high energy, it is only weakly 
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dependent on the region where the fit is performed. Alternate functional forms 
were also investigated, giving results consistent with those from Eq. 24. 



8.4-2 Error Correlations in Response Fit 

The error of the fit on the jet response measurement dominates the systematic 
uncertainty over large regions of jet energy and pseudorapidity. In this section, 
the point-to-point correlations in uncertainty due to the fit errors described 
in the previous section are derived. 

The error bands in Fig. 32 show the 68%-probability response at each energy 
based on fluctuations of all fit parameters. However, the parameter set that 
produces a given variation at one energy need not produce the same variation 
at any other energy in either magnitude or sign. In general, no single set 
of fit parameters may exist to produce the 68%-probability curves. But it is 
reasonable to expect that errors for two points of similar energy should be 
largely correlated. The correlations in the fit error are calculated as follows: 

• A grid of parameter sets (Ngrid) is generated to define the ^ xLin + 3-5 
volume. This volume contains all fit parameter set fluctuations correspond- 
ing to a 68% probability content. Each parameter set defines a response 
function contained within the bands shown in Fig. 32. 

• The parameter sets are used, noting the variation in response at 11 values 
of uncorrected jet energyQ between 10 and 500 GeV. The correlations in 
response for each pair of energy values are calculated from the parameter 
sets in the grid, which are used to define a correlation matrix for the 11 
energy values. Each element of the matrix {p{i,j)) is the standard correla- 
tion coefficient between the responses measured at each energy value, and 
defined as: 

p(.,.) = . , (26) 

ii:na'm,)-R(i)rz"R'(nM-mjmi 

where R{i) is the response for the i*^ energy bin. 

Table 5 shows the correlation matrix for the response fit to 0.7 cone jets. 
Correlations are illustrated graphically in Fig. 33, where four rows of the 
matrix are plotted, showing the error correlations relative to the errors at 20, 
50, 100, and 500 GeV respectively. 



In this context uncorrected energy means energy not corrected for response. The 
Iow-£^T bias, offset, and 77-dependent corrections have already been applied. 
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8.5 Shower Containment 



The D0 central and end calorimeters are more than 7.2 and 8 interaction 
lengths (Aint) thick, respectively. Although the EC is thick enough to con- 
tain all expected showers, there may be some leakage out of the CC for very 
high energy jets. Because the energy not contained within the calorimeter 
contributes to the measured response in principle corrects for leakage. 
The response of jets with energy less than ~ 100 GeV is measured using the 
CC and, therefore, includes the containment correction. For jets with energy 
^ 100 GeV, the response is determined from EC data. Measuring the high en- 
ergy response in the end cap calorimeter, with full containment, and applying 
it to the central calorimeter, where there could be some leakage, could cause 
a bias at high energy. The cryostat factor correction does not account for this 
bias, as it is measured from jets with energies ~ 100 GeV. The effect of shower 
containment on the response measurement at high energies is evaluated using 
a simulation based on Monte Carlo and experimental data. 

The NuTeV collaboration has measured the energy loss of charged pions as 
a function of number of interaction lengths in their calorimeter [19], which 
consists of stainless steel (absorber) and scintillator (active medium). For pions 
of various incident energies. Fig. 34 shows the fraction of energy deposited 
beyond a certain calorimeter depth as a function of that depth in units of 
interaction length, Aj^x- The energy loss as a function of Aint does not depend 
strongly on the type of absorber. The NuTeV data can, therefore, be used for 
this study despite the difference in the calorimeter compositions. 

The thickness of the D0 central calorimeter is determined from a GEANT 
simulation of the D0 detector. Figure 35 shows the depth of the CC as a 
function of t] in units of Aint- By definition, Aint depends on the proton 
cross section for inelastic collisions with nuclei of the absorption material. 
The uncertainty in this cross section accounts for the difference between the 
two sets of calorimeter thicknesses displayed in Fig. 35. GEISHA and fluka 
are different simulations of physics processes implemented in GEANT [9,20]. 

The energy containment of jets is modeled using particle level HERWIG jets 
for \ri\ < 0.7. The number of interaction lengths each particle traverses is de- 
termined from the 7] of each particle in the jet and the information displayed 
in Fig. 35. Figure 34 is used to determine the fraction of energy of each parti- 
cle that is contained in the calorimeter. All strongly interacting particles are 
treated as charged pions. Electromagnetically interacting particles (7r°,e,7) are 
fully contained. The energy contained within a jet is then compared to the 
total particle level jet energy to measure the fraction of energy escaping the 
calorimeter. 
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The results of the simulation arc shown in Fig. 36. The most conservative 
of the CC depth estimates (7.2 Ajnt at r]=0) is used. The data have been 
normalized to 1.0 at 100 GeV in order to evaluate only the energy loss with 
respect to 100 GeV jets. The effect on the measured response due to energy 
not contained by the CC is less than 0.5% for Et <, 400 GeV. 

The thickness of the calorimeter is known with a precision of ~ 0.4Aint- 
Different thicknesses between 7.2 and 7.6Aint are, therefore, modeled. While 
the absolute energy contained varies between the different models, the rel- 
ative change as a function of jet energy is constant. The simulation is also 
repeated using a parameterization of energy loss of charged pions as a func- 
tion of Aintj from Ref. [7]. The results are consistent with those obtained 
from the NuTeV parameterization. These results determine the assignment of 
a 0.5% uncertainty on the effect of shower containment on the CC response 
at high energies. 

8.6 Effects of Calorimeter Acceptance 

The event J^j' measured from energy deposited in the individual calorimeter 
cells. The sum is performed over all cells within the range |?7| < r/nm with r/um = 
4.5 {9 ^ 1°). To study the effects on i?jet of the large but finite calorimeter 
acceptance in the measurement, the recalculated as a function of 

77iiin in a 7-jet Monte Carlo sample (herwig-SHOWERLIB). From the measured 
Rjct versus ri^^^, it is possible to extract the effect of the limited acceptance by 
extrapolation to r/um — > 00. 

The calorimeter acceptance bias is negligible above i?jet ^ 50 GeV, given that 
-Rjct is independent of ?7iini above \ri\ = 4 (Fig. 37). Below this energy threshold, 
the bias is very small. A 0.5% error is assigned below Exjet = 15 GeV, which 
decreases linearly to zero at 50 GeV. 

8. 7 Summary of the Systematic Errors in Rjet 

In addition to the fit error, there is a 0.5% physics background error on Rj^t, 
as discussed in Section 6. The calorimeter acceptance error is 0.5% below 
15 GeV and decreases linearly to zero at 50 GeV. At very low E^s, below 
~ 20 GeV, the 1ow-£'t bias error quickly becomes the dominant uncertainty. 
Above Et — 100 GeV in the CC, the leakage uncertainty increases linearly 
from zero to 0.5% at 400 GeV. 

There are two biases present in the response measurement: the topology bias 
and the instrumental background bias. Due to the topology bias, the measured 
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response is 0.5%-l% low; this effect is constant as a function of jet energy (Fig. 
13). In addition, a constant bias of approximately 0.7% towards a higher re- 
sponse is observed due to instrumental background. It is difficult to accurately 
measure and correct for these biases. The two effects, however, approximately 
cancel each other, and therefore no correction is made. A 0.5% error is included 
to account for any residual bias. The leakage uncertainty is added linearly to 
this residual bias, before all the jet energy scale error components are added 
in quadrature. 

The i?jet correction is applied to ^/s — 630 GeV data as measured in the 
y/s = 1800 GeV 7-jet sample. While /?jct is expected to be the same at both 
ccntcr-of-mass energies, it is measured more accurately from the high statistics 
sample at y/s = 1800 GeV. Verification studies [4], based on independent 
measurements of Rj^t in both data sets, show that the response is the same 
to within 0.3%. This number is quoted as an additional uncertainty for jets 
produced during the low center-of-mass energy run. 



9 Showering Correction 

The last correction applied is the showering correction. It compensates for 
the net energy flow through the cone boundary during calorimeter showering. 
Ideally, the jet energy correction should scale the energy of the reconstructed 
jet to the particle level. As the particles comprising the jet strike the detector, 
they interact with the calorimeter material, producing a wide shower of addi- 
tional particles. Some particles produced inside the cone deposit a fraction of 
their energy outside the cone as the shower develops, and vice versa. 

It is not possible to determine this showering effect directly from calorime- 
ter data. Energy outside the cone may be associated with gluon emission or 
fragmentation outside the cone at the particle level ("physics out-of-cone" ) . 
Alternatively, such energy may be related to underlying event, noise, or pile- 
up. The correction is derived using the energy per Ar] x Aip in the vicinity of 
the jet centroid (energy density profile) obtained from both data and particle 
level HERWIG Monte Carlo. The physics out-of-cone energy (from particle level 
Monte Carlo) is subtracted from the total out-of-cone energy (from data) . 

9.1 Method 

The description that follows is based on 7^ = 0.7 jets with physics |?7| < 0.4. 
Annuli or sub-cones of increasing sizes (steps of 0.1) are defined around the 
jet centroid. -E'jet(r) is calculated as the sum over all cells contained within a 
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sub-cone of radius r, with r < 2. Contributions from the underlying event, 
uranium noise, and pile-up are subtracted. 



The total out-of-cone ratio, Ftot, is defined as the energy deposited inside a 
large cone taken as the hmit of the cluster divided by the energy deposited 
inside the algorithm cone. Ftot is measured as: 

_ i^jet(r < 1.0) 

- E,et(r < 0.7) • ^^^^ 

The assumption, based on the measured energy per Arj x A(p in the vicinity of 
the jet centroid, is that the cluster does not extend beyond r = 1.0 for central 
jets [3] . Showers produced by central and forward jets inside the calorimeter 
cover approximately the same area in real space. Pseudorapidity space, how- 
ever, shrinks towards the direction of the beam pipe. The cluster boundary, 
therefore, cannot be the same for all pseudorapidity bins. For |?7| < 0.4, the 
limit is chosen at r = 1.0. This value is increased towards the beam pipe up 
to r = 1.6 for 2.5 < \ri\ < 3. 

The factor Ftot includes both showering loss and physics out-of-cone. The 
latter, denoted as -Fphy, is obtained using the same procedure as for Ftot from 
a HERWIG particle level jet sample reconstructed with the 0.7 cone algorithm. 



The showering correction factor Fgho is defined as: 

" ,t(r <0.7) + r 
Fjet(r < 0.7) 



^s.o = ^44^^^^> (28) 



where Fgho is the amount of energy associated with particles emitted inside 
the cone at the particle level, but deposited outside the cone in an annulus 
defined by 0.7 < r < 1.0 at the calorimeter level. In the same way, Ftot and 
Fphy can be written as: 



^ Fjet(r < 0.7) + Fphy(r > 0.7) + F,ho 
£^jet(r < 0.7) 

p _ Fjet(r<0.7) + Fphy(r>0.7) 
£^jet(r < 0.7) 



(29) 
(30) 



where Fphy is the energy associated with the physics out-of-cone. Fgho can be 
expressed in terms of Ftot and Fphy as: 



sho 



i^tot - i^phy + 1 . (31) 
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In this paper, the fraction of the shower contained in the TZ—0.7 cone, S — 
l/Fgho, is measured. 

In the derivation of the showering correction described above, the offset is 
subtracted from the energy density profiles measured in data. After this sub- 
traction, the energy density is still not zero outside the cluster limits. This 
contribution to the energy density profiles is small and approximately a con- 
stant function of pseudorapidity. It probably comes from particles produced 
far from the jets (in a dijet events the two clusters are color connected) and 
constitutes one of the largest sources of systematic uncertainty at low energies. 
This baseline is subtracted from both the data and the particle level Monte 
Carlo, as the best solution to avoid the bias in S arising from differences in 
the baselines observed in both samples. 



9.2 Results 



Figures 38 and 39 show S versus £jet for the 0.7 cone over all physics pseudo- 
rapidities. The solid curves are fits to S versus E^ct- In the low energy range, 
the data are fit to either a logarithmic {a + b ■ InE'jct) or a linear (c + ■ -Ejct) 
function. At high energies, S is well described by a constant. The errors at 
low energy are dominated by the offset and baseline subtraction, and at high 
energy by fit errors due to poor statistics. The bands between dotted lines 
account for the total uncertainty. At high energies, the total uncertainty is 1% 
(4%) in the central (forward) regions for 0.7 cone jets. At very low energies, 
the errors increase up to 1% (10%) in the central (forward) regions for the 
same cone size. Figure 40 displays the parameterizations of S versus Ej^t for 
four cone sizes TZ —1.0, 0.7, 0.5, 0.3 and different rj regions. At high energies, 
the error for 0.3 cone jets is 2.5% (5%) in the central (forward) region. The 
uncertainty increases up to 2.5% (10%) in the central (forward) region at very 
low energies. Note that although the errors increase rapidly for very low energy 
jets in the forward region, they are not as large in the E^ range of interest. 
This is important because the significant variable in physics analyses is not 
energy but Et. 

Showering losses depend on the jet energy profiles in rj-(p space. Different 
profiles would translate into different responses. Given that the Rj^t measured 
from the ^/s — 1800 and 630 GeV samples are consistent, showering losses do 
not depend on the pp center-of-mass energy for jets of the same energy and 
pseudorapidity. 
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10 Monte Carlo Studies 



The Monte Carlo analysis described in this section serves two purposes. First, 
it provides the jet energy scale to correct Monte Carlo jets processed through 
the D0 shower library simulation (showerlib). Second, it proves that the 
method used to derive the correction achieves its purpose within errors. The 
analysis is based primarily on a sample of HERWIG direct photon events pro- 
cessed through SHOWERLIB for simulation of particle showers. 

10.1 Jet Energy Scale 

After application of the \ow-Et bias, offset, cryostat factor, and r^-dependent 
corrections, the Monte Carlo response in the CC is determined using the same 
procedure as for the data. Figure 41 shows R^^t versus E'jet for jets recon- 
structed with the 0.7 cone algorithm, along with the associated error band. 
The response is now uniform over the whole detector. As mentioned in Sec- 
tion 8, the shape of the response obtained from a HERWIG sample processed 
through SHOWERLIB is different from the response measured from data. In 
addition to a difference in overall normalization, the Monte Carlo response 
increases less rapidly, remaining nearly constant above 150 GeV. 

Some sources of error in the data analysis do not contribute to the Monte Carlo 
jet scale uncertainty. For example, uranium noise, pile-up, and multiple pp in- 
teractions are not modeled in the Monte Carlo sample. In addition, physics 
backgrounds in Monte Carlo are limited to diphoton events, which produce 
a negligible effect on the response measurement. Luminosity and multiple in- 
teraction cuts are not needed, because only single interactions are generated. 
There is also no instrumental background in the Monte Carlo. The topology 
bias is on the order of 1% at 30 GeV and becomes neghgible above ~ 200 GeV. 
Limits in the calorimeter acceptance contribute a small error to i?jet- Finally, 
because the full jet shower is contained within the calorimeter in the SHOW- 
ERLIB approximation, the shower containment error is negligible. 

10.2 Closure Tests 

The Monte Carlo sample provides an opportunity to verify the method used to 
derive the jet energy scale correction. This "closure" test directly compares the 

corrected jet energy with the energy of the associated particle jet. The ratio of 
these two quantities should be unity for all values of E^^^^ and pseudorapidity. 
For the closure test, the underlying event contribution is subtracted from 
both calorimeter and particle jets. Split and merged jets are removed from 
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the sample, because splitting and merging is not implemented in the particle 
level algorithms. 

Figures 42 and 43 show the ratio of calorimeter to particle jet energy before 
(open circles) and after (full circles) the jet scale correction, as a function 

of Ef^^^ and particle level rj. Statistical and systematic uncertainties are also 
shown. Within errors, the ratio is consistent with unity after all corrections. 
Thus, on average, the correction scales the energy of a 0.7 calorimeter jet to 
the energy of the associated particle jet to within ~0.5% for |?7| < 0.5. 

The closure test is also performed for other cones and pseudorapidity bins. 
Figure 44 shows the EJ^^'^yEf^^^ ratio versus £^Tjet calculated for the 0.7 cone 
jets in seven rj bins within the range \r)\ < 3. All ratios are consistent with 
unity within the total systematic uncertainties. 



11 Summciry and Conclusions 

The energy calibration was performed for jets observed in pp collisions with the 
D0 detector at Fermilab. The work described in this paper is based primarily 
on data taken by D0 during the 1992-1996 pp collider run. Test beam data 
and Monte Carlo samples arc also used in some cases. The corrections were 
derived for two center-of-mass energies of the pp system, y/s = 1800 GeV and 
^/s = 630 GeV, and are valid in the range Exjet > 8 GeV and |?7| < 3. The jet 
energy scale compensates for spectator interactions, uranium noise, response, 
and showering loss. The response correction is small, due to the hermeticity, 
hnearity, and good e/n ratio of the D0 calorimeters. 

Figures 45-52 show the magnitude of the total correction and uncertainty as a 
function of jet energy for pseudorapidities of 0, 1.2, and 2. The contribution of 
the individual sources of systematic error are also shown. The total uncertainty 
is the sum in quadrature of the individual components. In most of the plots 
in this paper, the fixed cone algorithm with TZ = 0.7 is used as an example, 
the data correspond to y/s = 1800 GeV collisions, and the luminosity is set 
to 5 X 10^° cm~^sec~^. 

The overall correction factor to jet energy in the central calorimeter is 1.160 
and 1.120 at 70 and 400 GeV, respectively. The total uncertainties at the 
same energies are 0.015 and 0.023. At lower energies, larger pseudorapidities, 
and smaller cone sizes, the corrections and errors increase. The procedure is 
verified with a Monte Carlo simulation which shows that the jet energy is 
corrected by this procedure to the particle level to within the quoted errors. 
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Table 1 

Quality cuts on photon candidates to reject background. 



Variable 


Er-y <15 GeV 


Et^ > 15 GeV 


dE/dx 


< 0.8 or > 1.5 


< 0.6 or > 1.5 




< 0.25 or > 0.75 


< 0.1 or > 0.9 


Ctrk 


> 3.0 


> 3.0 


EMF 


> 0.9 


> 0.96 


/iso 


< 0.5 


< 0.15 
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Table 2 

Means and RMS values for the residuals in //-dependent response derived from 7-jet 
data. Residual is defined as the difference between the response measured from data 
after the ?7-dependent corrections and the ideal response. 



T] Region 


Mean 


RMS 


< 0.5 


0.0 


0.008 


0.5 < \r]\ < 1.0 


0.004 


0.002 


1.0 < \r]\ < 1.5 


0.010 


0.005 


1.5 < \r]\ < 2.0 


0.004 


0.004 


2.0 < |?7| < 2.5 


0.0 


0.009 


2.5 < \r)\ < 3.0 


0.010 


0.010 
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Table 3 

Means and RMS values for the residuals in 77-dependent response derived from jet- 
jet data. Residual is defined as the difference between the response measured from 
data after the ?7-dependent corrections and the ideal response. 



r) Region 


Mean 


RMS 


\r)\ < 0.5 


0.0 


0.005 


0.5 < \r]\ < 1.0 


0.007 


0.005 


1.0 < \r]\ < 1.5 


0.001 


0.024 


1.5 < I77I < 2.0 


0.010 


0.022 
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Table 4 

Parameters of the fit to Rjet for different cone sizes. 



Jet Response Parameters 


Cone 


1.0 


0.7 


0.5 


0.3 


a 


0.6739 


0.6802 


0.6807 


0.6822 




±0.0485 


±0.0515 


±0.0506 


±0.0555 


b 


0.0433 


0.0422 


0.0429 


0.0439 




±0.0211 


±0.0236 


±0.0235 


±0.0263 


c 


-0.0013 


-0.0013 


-0.0014 


-0.0015 




±0.0025 


±0.0027 


±0.0027 


±0.0031 
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Table 5 

Correlation matrix for error band in jet response correction {TZ = 0.7). 



Correlations for Fit Error (TZ = 0.7) 



E (GeV) 


10 


20 


35 


50 


75 


100 


150 


200 


300 


400 


500 


10 


1.00 


0.98 


0.71 


-0.37 


-0.70 


-0.65 


-0.39 


-0.13 


0.21 


0.38 


0.48 


20 


0.98 


1.00 


0.81 


-0.24 


-0.65 


-0.64 


-0.43 


-0.20 


0.11 


0.28 


0.37 


35 


0.71 


0.81 


1.00 


0.35 


-0.19 


-0.28 


-0.25 


-0.15 


0.00 


0.10 


0.15 


50 


-0.37 


-0.24 


0.35 


1.00 


0.83 


0.70 


0.50 


0.33 


0.09 


-0.03 


-0.10 


75 


-0.70 


-0.65 


-0.19 


0.83 


1.00 


0.97 


0.82 


0.62 


0.32 


0.14 


0.04 


100 


-0.65 


-0.64 


-0.28 


0.70 


0.97 


1.00 


0.93 


0.78 


0.51 


0.34 


0.24 


150 


-0.39 


-0.43 


-0.25 


0.50 


0.82 


0.93 


1.00 


0.96 


0.79 


0.67 


0.58 


200 


-0.13 


-0.20 


-0.15 


0.33 


0.62 


0.78 


0.96 


1.00 


0.94 


0.85 


0.79 


300 


0.21 


0.11 


0.00 


0.09 


0.32 


0.51 


0.79 


0.94 


1.00 


0.98 


0.95 


400 


0.38 


0.28 


0.10 


-0.03 


0.14 


0.34 


0.67 


0.85 


0.98 


1.00 


0.99 


500 


0.48 


0.37 


0.15 


-0.10 


0.04 


0.24 


0.58 


0.79 


0.95 


0.99 


1.00 
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Fig. 1. General view of D0 calorimeters. 
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Fig. 2. Diagram of a D0 calorimeter unit cell. 
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Fig. 3. Side view of the D0 calorimeters (one quadrant). The numbers indicate 
position in units of pseudorapidity. 
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Fig. 4. The ratio c/tt as measured from test beam data and determined from 
a Monte Carlo simulation. Note that the Monte Carlo e/vr ratio (open squares) 
changes faster and flattens out earlier than the test beam ratio (full circles). 
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Fig. 5. Measured test beam energy versus particle momentum for electrons and 
pions. The solid and dashed lines are fits to the electron and pion data, respectively. 
Good linearity between reconstructed and test beam energy is achieved with the 
D0 calorimeters. 
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Fig. 6. {Et\ — Et2)/{Eti + Et2) distribution from dijet events observed in the 
D0 central calorimeter. The data are well described by a Gaussian, showing the 
hermeticity and linearity of the calorimeters. 
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Fig. 7. Physics underlying event Et density, Due (Vs = 1800 and 630 GeV). The 
larger samples available for the low center-of-mass energy measurement explain the 
smaller statistical errors. 
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Fig. 8. Transverse energy density, Dq, contributed by uranium noise, pile- up, and 
additional hard interactions, for different luminosities in units of lO^'^ cm~^sec~^. 
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Fig. 10. Parametric simulation of the Rjet measurement. A fit to Rjet "measured" 
from the simulated data, when binned in terms of E' to remove resolution effects, 
agrees very well with the input response (top). The agreement between the fit and 
the input function is also excellent given different reasonable assumptions for the 
Et dependence of the 7 cross section, such as E^^, E^^, and E^"^ (bottom). 
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Fig. 11. Change in the measured i?jet as the cut is tightened, for three different 
E' bins. i?jet (std cuts) is the measured response from the 7-jet sample selected with 
the cuts listed in Tabic 1. An additional cut on (matrix test) is applied for the 
other bins to study the effect of highly EM jets on i^jet- 
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Fig. 12. i?jet as a function of the Aip cut threshold for 50 < .E' < 60 GeV. The 
circles correspond to an integrated distribution (lower Acp bins include the events 
of the higher A(p bins), while the open squares represent the differential distribution. 
The solid and dashed lines are parameterizations of the differential and integrated 
distributions. The cut is at A(p = 2.8 radians (arrow). 
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Fig. 13. Difference between i?jet at Aip = tt and A(p > 2.8 radians versus E'. The 
circles are the result from tfie integrated distribution and the squares from the 
differential distribution. 
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Fig. 14. i?jet versus E' in the CC and EC from a low luminosity sample 
(£ < 5 X 10^° cm~^sec~^). If the single interaction requirement is not enforced, 
the response distribTitions are biased towards lower values, due to the effect of 
misvertexing (full circles). The bias vanishes if only single interaction events are 
accepted (open triangles). 
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Fig. 15. Response as a function of luminosity for CC jets with E' = 35, 60, and 
120 GeV. The variable delta(0,5) is the difference in R^et measured from samples 
taken at luminosities of zero and 5 x lO^'^ cm~^sec~^. The single interaction require- 
ment is always enforced. 
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Fig. 16. Response as a function of luminosity for EC jets with E' = 120, 190, and 
290 GeV. The variable delta(0,5) is the difference in i?jet measured from samples 
taken at luminosities of zero and 5 x 10^°cm""^sec~^. The single interaction require- 
ment is always enforced. 
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Fig. 17. Rjet versus E' for a sample with no ??-dependent corrections. CC jets cor- 
respond to |r/| < 0.7, IC jets to 0.7 < |r/| < 1.8, and EC jets to 1.8 < |r/| < 2.5. 
The lowest E' point in the CC is affected by 1ow-^t resolution bias, as explained 
in Section 8. 
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Fig. 18. Cryostat factor as a function of E'. The solid line is a fit to a constant 
function. 
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Fig. 19. Response versus detector r] for jet-jet data (?7 > 0). The line is the fit to 
the ideal rj dependence. 
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Fig. 20. Response versus detector r] for 7-jet data. The dashed Hne is the fit to the 
ideal 77 dependence. 
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Fig. 21. Correction factor Fjj as a function of central jet Et for positive rj between 
0.8 and 1.2. The circles correspond to the 7-jet data and the triangles to the jet-jet 
data. The line is a linear fit through all the data points. 
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Fig. 22. Response versus 77 for 7-jet events after all corrections have been applied 
(including the energy dependent correction). 
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Fig. 23. Response as a function of from jet-jet data after all corrections have been 
applied (averaged over positive and negative 77). The line is the ideal 77 dependence. 
A resolution bias is evident at large 77 (outside the IC region), where no correction 
needs to be derived. 
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Fig. 24. Low-Et bias, i?bias) as a function of £^rjet (7^=0.7). The full circles denote 
the simulation and the open circles the data measurement. The solid and dashed 
lines are fits to the simulation and the data, respectively. The dotted band is the 
total systematic uncertainty in the simulation. 
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Fig. 25. Rjet versus E' after r/-dependent corrections are applied. CC jets correspond 
to \rj\ < 0.7, IC jets to 0.7 < \r]\ < 1.8, and EC jets to 1.8 < \rj\ < 2.5. 
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Fig. 26. -Ej.jgj'^ (after offset and ry-dependent corrections) versus E' for the 0.7 cone 
jet algorithm. The mismatch between CC and EC jets is due to showering effects. 
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Fig. 27. Comparison between i?jet measured from data and i?jet determined from a 
SHOWERLIB Monte Carlo sample iJZ = 0.7). The Monte Carlo response flattens out 
more rapidly and is nearly constant above 150 GeV. 
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Fig. 28. (a-c) Single particle responses for electrons and pions from test beam data, 
(d) Single particle response from Monte Carlo. The test beam responses in (a-c) are 
identical above 2 GeV; they differ in the assumptions made ioi E < 2 GeV where 
no data are available. 
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Fig. 29. Jet energy response determined using the projection fraction method 
from SHOWERLIB Monte Carlo (circles, triangles and stars) compared to the response 
obtained with the single particle convolution approach (dashed line). The solid line 
shows a fit to the points. 
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Fig. 30. The Monte Carlo response using the convolution method and the best test 
beam model for single particle response is shown with the stars. It is normalized to 
i?jet as measured from collider data, shown as open circles in the CC and full circles 
in the EC. 
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Fig. 31. The solid curves show the extremes of the three test beam response models 
normalized to the collider data. The star indicates an estimated response from the 
convolution method for 500 GcV jets. The inner error bar is based on the difference 
between the two extreme models. The full error bar includes all errors, as explained 
in the text. 
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Fig. 32. ii!jet versus energy on a logarithmic (top) and linear (bottom) energy scale 
for the 0.7 cone jet algorithm. The outer band shows limits on the measured response 
for jets based on the region in parameter space defined by the = X^in+^-S surface. 
This region corresponds to the 68% confidence region of parameter fluctuations from 
the nominal values. 
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Fig. 33. Error correlations for 0.7 cone response fit. Error correlation points are 
shown in four slices from the full correlation matrix. The lines connect the points. 
The four curves show the point-to-point correlation of fit errors relative to energy 
values of 20, 50, 100, and 500 GeV respectively. 
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Fig. 34. Fraction of pion energy escaping from the NuTeV calorimeter as a function 
of depth in units of interaction lengths. 
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Fig. 35. Depth of the D0 central calorimeter in units of interaction lengths. 
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Fig. 36. Fraction of the jet energy contained within the D0 central calorimeter 
{\r]\ < 0.7) as a function of jet energy. The data are normalized to unity at 100 GeV. 
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Fig. 37. i?jet in the CC as a function of the r]\\^ used in the Ifirj, calculation. The 
calorimeter acceptance bias is negligible above E'jct ^ 50 GeV, given that iijet is 
independent of the acceptance pseudorapidity threshold above \ri\ = 4. 
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Fig. 38. S versus Ejet (T^ = 0.7). Ejet is corrected by low-Ex bias, offset, and 
response. Tlie solid line is a fit to tfie data. The dotted lines sfiow tfie systematic 
uncertainty. 
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Fig. 39. S versus £jet = 0.7). E^^t is corrected by low-£^r bias, offset and response. 
The solid line is a fit to the data. The dotted lines show the systematic uncertainty. 
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Fig. 40. S versus £jet for different cone sizes in different ij regions. Only the param- 
eterizations to the data points are shown. 
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Fig. 41. ii!jet values for TZ = 0.7 Monte Carlo jets (full dots are CC jets, triangles 
EC jets and stars IC jets) along with the result of the nominal fit and the error 
band (herwig-showerlib). 



83 



^ 1.1 



LJ 



S 1.05 

E 



0.95 



0.9 



0.85 



0.8 



0.75 



<^ 6 ^ 



50 



100 



■f-f 



o o 



l?7l<0.5 0.7 cone 



• Corrected E" 
o Uncorrected E" 



jet 



jet 



150 



200 



250 300 350 400 450 

E^'jetCGeV) 



Fig. 42. Ratio of calorimeter and particle jet energy before (open circles) and after 
(full circles) the jet scale correction is applied (herwig-SHOWERLIb). The ratio is 



plotted as a function of E?^^^ for TZ = 0.7 Monte Carlo jets. Inner bars are statistical 
errors, and outer bars contain statistical and systematic errors added in quadrature. 
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Fig. 43. Ratio of calorimeter and particle jet energy before (open circles) and after 
(full circles) the jet scale correction is applied (herwig-SHOWERLIb). The ratio is 
plotted as a function of particle level r/ for 7?. = 0.7 Monte Carlo jets. A Et > 15 GeV 
cut is applied to remove the region where showering losses and Iow-Et bias effects 
dominate. Inner bars are statistical errors, and outer bars are statistical and sys- 
tematic errors added in quadrature. 
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Fig. 44. Ratio of calorimeter and particle jet energy after the jet scale correction 
is applied (herwig-SHOWERLIb). The ratio is plotted as a function of E^^^^^ for 
TZ = 0.7 Monte Carlo jets. The bars are statistical and systematic errors added in 
quadrature. 
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Fig. 45. Corrections and errors for r] = 0, TZ = 0.7. Top: Nominal, high (nominal 
+ (t), and low (nominal — a) correction factors. Bottom: high and low fractional 
errors. 
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Fig. 46. Individual components of the fractional jet scale error versus uncorrected jet 
Et for 77 = (7?. = 0.7): physics underlying event, uranium noise (zero suppression), 
and multiple interactions (Section 4); cryostat factor and ij dependence in the IC 
region (Section 7); low-E'^^ bias, response fit, physics background, and residual bias 
from event topology, instrumental background and shower containment (Section 8); 
showering (Section 9). Some components are zero because they do not contribute 
in that particular 77 region. 
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Fig. 47. Total Corrections and errors for |?7| = 1.2 (TZ = 0.7). Top: Nominal, high 
(nominal +cr), and low (nominal — cr) correction factors. Bottom: high and low 
fractional errors. 
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Fig. 48. Individual components of the fractional jet scale error versus uncorrected 
jet Et for rj = 1.2 (TZ = 0.7): physics underlying event, uranium noise (zero sup- 
pression), and multiple interactions (Section 4); cryostat factor and rj dependence in 
the IC region (Section 7); 1ow-£'t bias, response fit, physics background, and resid- 
ual bias from event topology, instrumental background and shower containment 
(Section 8); showering (Section 9). 
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Fig. 49. Corrections and errors for r] = 2 {TZ = 0.7). Top: Nominal, high (nominal 
+ and low (nominal — a) correction factors. Bottom: high and low fractional 
errors. 
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Fig. 50. Individual components of the fraetional jet scale error versus uncorrected jet 
Et ioi t] = 2 (TZ = 0.7): physics underlying event, uranium noise (zero suppression), 
and multiple interactions (Section 4); cryostat factor and 77 dependence in the IC 
region (Section 7); low-Ex bias, response fit, physics background, and residual bias 
from event topology, instrumental background and shower containment (Section 8); 
showering (Section 9). 
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Fig. 51. Corrections and fractional errors versus rj (jet Et = 20 GeV and TZ = 0.7). 
The total correction and error are both shown as well as the rj dependence of several 
individual components of the jet scale error: cryostat factor and ij dependence in the 
IC region (Section 7); residual bias from event topology, instrumental background 
and shower containment (Section 8); showering (Section 9). 
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Fig. 52. Corrections and errors versus rj (jet Et = 100 GeV and TZ = 0.7). The total 
correction and fractional error are both shown as well as the r] dependence of several 
individual components of the jet scale error: cryostat factor and r] dependence in the 
IC region (Section 7); residual bias from event topology, instrumental background 
and shower containment (Section 8); showering (Section 9). 
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